Fluctuation Pressure of a Membrane Between Walls Through Five Loops 
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An earlier four-loop calculation of the fluctuation pressure of a fluid membrane between two infinite 
walls is extended to five loops. Variational perturbation theory is used to extract the hard-wall limit 
. from perturbative results obtained with a smooth potential. Comparison with a structurally similar 

£f~) ' quantum mechanics problem of a particle in a box is used for an alternative way of extracting the 

membrane pressure and also to estimate the quality of the results. Our values lie above the best 
available Monte Carlo data. 

O 1 PACS numbers: 05.40.-a, 46.70.Hg, 87.16.Dg, 05.10.-a 

I. INTRODUCTION 

53- 
(— > , 

c/3 i The dominant repulsive force between layered chemical and biological systems, called membranes, is given by 
thermal out-of- plane fluctuations [f , 2]. In the absence of tension, these membranes are called fluid and the fluctuations 
are controlled by the membranes' bending rigidity. There have been various theoretical approaches to compute the 
pressure of a single membrane between walls [1, 3-7] or of a stack of membranes [1, 3-5, 8]. These situations are also 
' interesting statistical mechanics problems. 

Here we are concerned with the pressure generated by the bending fluctuations of a fluid membrane between two 
Q ■ infinitely extended parallel walls, which has the form [1] 

o : 

>■..... 

1 where k is the bending rigidity of the membrane, d is the distance between the walls and a is a factor that we wish 
f*"*. ' to compute. Estimates of a have been ranging over the years from theoretical approximations a m 0.0242 by Hclfrich 
[1] and a w 0.0625 by Janke and Kleinert [3] through Monte Carlo estimates a — 0.079 ± 0.002 by Janke, Kleinert 
and Meinhart [4] and a — 0.0798 ± 0.0003 by Gompper and Kroll [5] and a recent theoretical estimate a w 0.0797 by 
Bachmann, Kleinert and Pelster [7]. 

In [7], the result was obtained by replacing the hard walls by a tan 2 potential, whose prefactor was sent to zero at the 
end of the calculation to recover hard walls. This corresponds to the strong-coupling limit a = lim^oo a(g) of a loop 
expansion of a(g), where 1/g 2 is proportional to the prefactor of the potential. To achieve the necessary resummation, 
variational perturbation theory (VPT) [9] was used. This technique has been successful also in other situations where 
the strong-coupling limit of an asymptotic weak-coupling series is sought, e.g. when computing critical exponents from 
(j) 4 field theory models [10]. In this work, we extend the four-loop calculation of [7] to five loops. 

Our work is structured as follows. In sec. II, we model the hard walls with two different potentials and give 
the perturbative results in sec. III. In sec. IV, we follow [7] and use VPT to estimate the strong-coupling limit 
corresponding to hard walls. In sec. V, we consider the two different potentials to model the walls for a quantum 
mechanical (QM) particle in a box. The problem of finding the ground state energy for the QM problem is identical 
to computing the partition sum of a string between walls modeled by the same potential [6, 11, 12]. This problem 
in turn is structurally equivalent to finding a(g) in the membrane problem. Only for one of these potentials, the 
solution is known exactly [6, 12]. Although the potentials are very similar in the region of interest, their behavior 
under resummation with VPT is rather different. This will be used to judge the quality of the results of using VPT 
for the membrane problem. In sec. VI we force a(g) for the membrane problem to be identical to a(g) in the solvable 
QM problem by choosing the potential appropriately and extract a from determining the potential's singularities. In 
sec. VII we discuss our results. 
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II. MODELING OF THE BOUNDARY CONDITIONS 

Consider a tensionless membrane between two large flat parallel walls of area A separated by a distance d. In the 
harmonic approximation, which we are considering throughout, the curvature energy is given by 



E = \ l A d 2 x[d\{x)]\ (2) 



where n is the membrane's bending rigidity and ip is a field that describes the membrane's position between the walls, 
which are located at ±d/2. The ci-dependent part fa of the free energy density of the system at temperature T is 
given by the path integral 



ex K-^) = n/_r #(:c)exp (-^ 



(3) 



The pressure is then obtained as 



P= -^d (4) 



and has the form (1) [1, 3] and our task is to find the constant a. 

The difficulty in computing the path integral (3) consists in implementing the restriction — d/2 < ip < d/2. We 
follow [6, 7] and add a potential term m 4 d 2 J d 2 xV((p/d) to E, where V has a sufficiently strong singularity at ±1/2, 
expand the potential V in a Taylor series in ip and drop the restriction on <p. At the end of the calculation we take 
the limit m — > 0. We consider the potentials 



2-k 1 co^yjtz) 



and 



(l + 2z) 2 (l-2z) 2 



(6) 



The potentials have in common that they have quadratic divergences at ±1/2 and that their quadratic term in a Taylor 
expansion is normalized to z 2 /2. V c is related to the potential V t = (2tt 2 )^ 1 tan 2 (7ra;) used in [7] by V c = (2tt 2 )^ 1 + V t . 
For the resummation procedure employed in [7], V c and V t yield identical results and we will therefore recover the 
four-loop result reported there. For the other procedures used here, V c is better suited than V t . 

Since the functional form of p in terms of n, d and T is known and since we are going to differentiate only with 
respect to d, we will set UbT = n = 1 in the sequel. The energy functional may then be expanded as 



E 



jd 2 x\± [d 2 p(x)] 2 + \m\{ X f + m 4 e d 2 + m 4 £ e 2k d 2 ^ p(x) 2k j , (7) 



where the e 2 fc are the expansion coefficients of the potential. 

The path integral can now be evaluated in a loop expansion [6, 7]. The resulting Feynman diagrams, including 
their combinatorial factors, are obtained from recursion relations, whose derivation is delegated to appendix A. The 
evaluation of the associated momentum integrals is detailed in appendix C. 

III. PERTURBATION THEORY 

The diagrams labeled L—n (nth L-loop diagram) of appendix A correspond apart from combinatorial factors 
cl-u and coupling constant factors gL~n to integrals in a;-space. Since the diagrams are connected and because of 
translational symmetry in the infinite- wall limit, we may split off a factor A from each diagram and represent the 
remainder in momentum space. A line represents then a propagator 

A(p 2 , m 2 ) = 4 ^ 4 - -lj ( 2 ' 1 . 2 - 2 -. 2 ) , (8) 
+ to 4 2m z \p z + im z p z — im A J 
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while the integration measure over all independent momenta is 

d 2 p 



(2^ 



with momentum conservation at each vertex. A vertex with 2k legs represents a factor 

-m 4 d 2 ^- k h 2k . 



(9) 



(10) 



The sum of all diagrams corresponds to the negative of the free energy density f m , where the index refers to the 
presence of a non-zero m and lim m ^o fm = fd- 

In the sequel, a diagram represents only the corresponding momentum space integral which we call Il-ui i- e - we 
split off not only a factor A, but also the combinatorial factor cl_„ and the — £2fc-part of the factors (10), which we 
collect into gL-n- Then f m has L-loop expansions 



1 L 

■> m = ^2 1^ aig ' 



1=0 



with 



and a coupling constant 



— - ^ gL-nCL-nlh-n 



m 2 d 2 



(11) 



(12) 



(13) 



In Table III, we give gL-n, CL-n an d Il-u through five loops. For instance, the resulting zero-, one- and two-loop 
contributions are 



1 

o 

Through five loops, we get for the potentials under consideration 
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(14) 



L 


aL for V c 


a L for V a 





0.0506606 


0.0416667 


1 


0.125000 


0.125000 


2 


0.154213 


0.156250 


3 


0.105998 


0.102307 


4 


0.026569 


0.028101 


5 


-0.034229 


-0.031426 



(15) 



IV. a FROM STRONG-COUPLING VARIATIONAL PERTURBATION THEORY 

The (i-dependent part of the free energy has for m 2 = the form / = Aa/d 2 , where the factor 4 ensures consistency 
with (1). Our task is to find an approximation to the strong-coupling limit a = lim^oo a(g) with L-loop expansions 
of a(g) given by 

1 L 

a (9) = 73 ^2 a i9 L , (16) 

g i=o 

with the knowledge of only the first few a;. We will assume that a(g) has a strong-coupling expansion 



= I E a 'm9- 2m/q 
m=0 



(17) 



4 



with an additional parameter q. Then the problem has the following form: Given a function /(<?) = 4a(g) with L-loop 
weak-coupling expansions 



fL(g) = 9 r J2fl"9 l 



(18) 



and assuming strong-coupling expansions 



M 



f M (g) = 9 P/9 X fm9~ 2m/q , 

m=0 



(19) 



we are interested in finding /J, p and q. Assuming a thermodynamic limit for the problem at hand means setting 
p = 0. Then a exists and is given by a = /d/4. In [7] it was additionally assumed that q = l, which is motivated by 
a similar QM problem, see sec. V. 
In VPT [9], we replace in (18) 



g l+r ^(tg) l+r l(^) 2/q + t 



iJ /q (tgy^Si + t 



2/q 



2/9 



[p-(l+r)q]/2 



1 ) [p-(l+r)q]/2 



1 



(20) 



reexpand the resulting expression in t through t L+r , set t = 1 and then optimize the resulting expression in g, where 
optimizing refers to the principle of minimal sensitivity [13] and in practice means finding appropriate stationary or 
turning points. That is, we replace 



J+r 



L-l 



g) 



k=0 



2/9 



-i k 



(21) 



and optimize the resulting expression in g. For V c and V a , where r = —2, we obtain with p = 



1 

4^3 



L-l 



.1=0 



k=0 



(2-09/2 
k 



(22) 



the 



as the L-loop variational approximation to a. This expression also holds for V t , where r = — 1. For q = 1, 
expression in square brackets is independent of ao, which is why we reproduce below the results of [7]. 

If we do not want to make assumptions about q for f m , we can determine it self-consistently by first treating 
din f^/ ding in VPT, since it has the same q as f m and since 



lim 



ding 



p 
<1 



(23) 



with p = by assumption of a thermodynamic limit. That is, we resum the expansion of din f^/ ding as detailed 
above and tune q such that optimization with respect to g leads to din f^/ ding = 0. 

A similar QM problem (see sec. V below) leads us to try q = 1. Let us consider the different potentials for modeling 
the walls that enclose the membrane with this assumption. The loop orders 0-2 do not admit a variational solution 
and we therefore take the perturbative results as our best approximation. Then the loop orders and 1 yield zero for 
a, since they contain only negative powers of g, and the two-loop result is a 2 = a 2 /4. The results through five loops 
are 



L 


a for V c 


a for V a 


2 


0.038553 


0.039063 


3 


0.073797 


0.073688 


4 


0.079473 


0.079422 


•5 


0.081354 


0.081345 



(24) 
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with the results for V c through four loops coinciding with those reported in [7]. An extrapolation of the results (24) 
suggests a value of a between 0.0820 and 0.0825. 

The results from determining q self-consistently as described above are through five loops 





V c 


v a 


L 


q 


a 


q 


a 


3 


0.38124 


0.093076 






4 


0.56789 


0.095830 


0.46463 


0.098222 


5 


0.73907 


0.090983 


0.74209 


0.090321 



The values are compatible with convergence towards q = 1 and with the a- values for q = 1, but convergence is too 
slow for any quantitative use. 



V. QUANTUM MECHANICAL PARTICLE IN A BOX 

A one-dimensional problem similar to the two-dimensional case above is finding the ground state energy of a QM 
particle in a one-dimensional box [12]. The Euclidean path integral to be computed is 



r+d/2 

exp(-T£ (0) ) =TT / dip(t)exp(-E) (26) 

\ J-d/2 



with 



1 



-T/2 



E=^l dt<p(tf, (27) 

L J -T/2 

where T is the total interaction time, being equivalent to the area A in the membrane case. In the large-T limit, 
is the ground state energy, for which we will test our approximation methods. Its exact value is 



E<°> - ij. ,28) 



Again, we model the walls with a potential, 



f+T/2 

E= dt 

J -T/2 



l -^{tf + m 2 d?V^{t)/d) 



(29) 



and as in the membrane case for f m /(kBT), we can write down a loop expansion for E(°\ After modifying the 
Feynman rules according to 

/IW3 - 4 ^-"«— -» v<i -' ) - < 3 »> 



p 4 + TO 4 p 2 + TO 2 ' 

and defining a (g) and g for the QM problem by 



not only can a(g) be expanded as in (16), but due to the simple relation (B5) between one-loop integrals in the 
membrane and QM cases, all diagrams that separate into one-loop integrals give the same contribution to a(g) in 
both cases [6]. It follows that for any given potential, ao, ai and 02, which involve at most one-loop topologies, are 
identical in the QM and the membrane problem. 

For V c , the exact ground state energy is known for any m and d [12], 



2d 2 \ ttV 2 7r 2 5 V 64 ) 2d 2 \2 w 4 g 2 2\ ir 4 g 2 



m 7T 2 / 16 1 4 / TTV \ 7T 2 / 1 16 1 / 6T , 
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The limiting value for g — > oo is in each case 

„2 



a = ^—~ 0.0771063. (33) 
128 



The coefficients a; and a' m in the weak- and strong-coupling expansions (16) and (17), respectively, can be obtained 
to arbitrary order simply by Taylor-expanding (32). Note how this implies q = 1 in VPT. This is the reason why we 
used q = 1 in VPT for the membrane problem. 

While for general potentials, the ground state energy cannot be computed exactly, it is possible to compute all 
Feynman diagrams analytically (see appendix B). Alternatively, it is possible to compute the coefficients ah to arbi- 
trary order by generalizing [12] the Bender- Wu recursion relation for the anharmonic oscillator [14]. The generalized 
relation reads (correcting some typos in [12]) 

n n—1 

4jc nj = 2(j + l)(2j + l)c nJ+1 - ^2(-l) k e 2 k+2C n -k,j-k-i - 2 ^ c k ic n - k ,j, l<j< 2n, (34) 

fc=i fe=i 

cqo = 1 and in all other cases c n j = 0. The are then given by 

a L = (-^j cl-1,1 L>2. (35) 

The results of carrying out VPT through 20 loops for the potentials V c and V a are collected in Table I and illustrated 
in Figs. 1 and 2. For fixed q = 1 we get exponentially fast convergence towards the exact value of a for V c . For V a 
no convergence is obvious, although the values obtained through the order considered are not far from the exact a. 
Essentially the same is true when determining q self-consistently, except that the convergence towards the exact value 
of a is delayed as compared to taking q = 1. For V c , q = 1 is approached exponentially fast, while for V a , q > 1 seems 
preferred at higher orders. 




FIG. 1: Quantum mechanical particle in a box. a as a function of the loop order L for q — 1 for V c (dashed line), V a (solid 
line) and the exact result (horizontal line). Note how the convergence towards the exact result is exponentially fast for V t , 
while questionable for V a - 

It is likely that the inferior convergence behavior for the potential V a originates in our missing understanding of the 
analytical structure of as a function of g. It is possible that the strong-coupling behavior is not of the form (19) 
or that the strong-coupling expansion has a zero radius of convergence. Numerically, the deviations of the coefficients 
aL from those for V c are relatively small in low orders, in particular the deviation from = for small even L > 2. 
Note how this is very similar to the results of the membrane loop expansion (15). Another likely similarity between 
the membrane problem and the QM problem with potential V a is the factorial growth of with L for large L. In 
this respect, using V c in the QM problem is very special, as already noted in [12], and it appears likely that the a_L 
grow factorially for the membrane problem for both V c and V a - Note, however, that the results for a with q = 1 and 
also for a with self-consistently determined q improve with increasing L as long as the a^ do not significantly grow. 
We will come back to this point in sec. VII. 
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V c 


Va 






q = l 


q self. -cons. 




g = l 


q self.-cons. 


L 


ah 


a 


q 


a 


aL 


a 


q 


a 





0.0506606 








0.0416667 








1 


0.125 








0.125 








2 


0.154213 


0.0385530 






0.15625 


0.0390625 


0.309401 




3 


0.0951261 


0.0719411 


0.605551 


0.0836038 


0.0911458 


0.0717445 


0.630222 


0.0819600 


4 





0.0758821 


0.850234 


0.0807166 


0.00325521 


0.0758318 


0.805894 


0.0816667 


5 


-0.0361959 


0.0767518 


0.931591 


0.0787187 


-0.0340667 


0.0767990 


0.920850 


0.0789522 


6 





0.0769910 


0.966170 


0.0778393 


-0.012597 


0.0770078 


0.975808 


0.0775787 


7 


0.0275454 


0.0770659 


0.982590 


0.0774492 


0.0421369 


0.0770326 


0.994979 


0.0771334 


8 





0.0770913 


0.990852 


0.0772701 


0.0400356 


0.0770777 


0.975795 


0.0774386 


9 


-0.0262028 


0.0771005 


0.995143 


0.0771857 


-0.164914 


0.0771337 


0.957841 


0.0778059 


10 





0.0771040 


0.997410 


0.0771451 


-0.207989 


0.0771252 


1.000970 


0.0771197 


11 


0.0279168 


0.0771054 


0.998617 


0.0771254 


1.56427 


0.0770648 


1.018330 


0.0767930 


12 





0.0771059 


0.999262 


0.0771157 


2.21468 


0.0770447 


1.008300 


0.0769460 


13 


-0.0318674 


0.0771061 


0.999607 


0.0771110 


-25.1291 


0.0771620 


1.005830 


0.0770638 


14 





0.0771062 


0.999792 


0.0771086 


-43.0543 


0.0771831 


0.979132 


0.0775060 


15 


0.0381093 


0.0771063 


0.999890 


0.0771074 


585.908 


0.0771362 


1.028460 


0.0766838 


16 





0.0771063 


0.999942 


0.0771069 


1288.21 


0.0771241 


1.038400 


0.0765030 


17 


-0.0471274 


0.0771063 


0.999969 


0.0771066 


-18478.5 


0.0771695 


1.029930 


0.0766691 


18 





0.0771063 


0.999984 


0.0771064 


-53154.3 


0.0772534 


1.037170 


0.0765617 


19 


0.0597739 


0.0771063 


0.999992 


0.0771064 


753376 


0.0772069 


1.024240 


0.0768319 


20 





0.0771063 


0.999996 


0.0771063 


2833593 


0.0772520 


1.051640 


0.0762912 



TABLE I: Determination of a (exact value a = 0.0771063 . . .) for QM particle in a box through 20 loops. Values of a for both 
potentials with q — 1 and values of q and a for both potentials with self-consistently determined q are shown. 



VI. MEMBRANE PROBLEM WITH a(g) FROM PARTICLE IN A BOX 

Let us compare the values for the QM expansion coefficients and the corresponding membrane coefficients using 
the potential V c : 



L 


QM 


^memb 
a L 


„mcmb „QM 
a L — a L 


1 


0.125 


0.125 





2 


0.154213 


0.154213 





3 


0.095126 


0.105998 


0.010872 


4 





0.026569 


0.026569 


5 


-0.036196 


-0.034229 


0.001967 



(36) 



We see that the relative difference through the order considered is small when both a™ e and a£ are nonzero. 
This is the motivation to carry out a different procedure for finding a from the loop expansion. Instead of asking 
directly what a is for a given potential for the membrane case, we slightly modify the order by order such that the 
expansion of a(g) is identical to that of the QM case with potential V c and ask where the resulting potential has the 
nearest singularity. The scaling relation / oc 1/d 2 when m 2 — allows us then to recover a for the membrane case. 
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FIG. 2: Quantum mechanical particle in a box. (a) Self-consistently determined g as a function of the loop order L from V c 
(dashed line), V a (solid line) and q = 1. Note how the convergence towards q = 1 is exponentially fast for V c , while for V a , no 
convergence is obvious, although q is around 1. (b) a as a function of the loop order with self-consistently determined q for 
V c (dashed line), V a (solid line) and the exact result (horizontal line). Note how the convergence towards the exact result is 
exponentially fast for V c , while questionable for V a . 



The expansion coefficients of the potentials are 





QM: V c 


memb. 




0.0506606 


0.0506606 


£2 


0.5 


0.5 


£4 


3.28987 


3.28987 


£6 


18.3995 


18.0284 


es 


94.6129 


89.5702 


£10 


462.545 


419.568 



(37) 



Let us instead investigate the expansion of 

oo 

l/v^j^V 21 , (38) 

k=0 

since for the QM case, we have 1/ y / 2n 2 V(x) = cos(nx). We can expect a good approximation for the location of the 
singularity of V if this singularity is of the quadratic type as in V c and V a . The expansion coefficients of the quantity 
(38) for the QM and membrane cases are 





QM 


memb. 




1 


1 


v 2 


-4.93480 


-4.93480 


v 4 


4.05871 


4.05871 


v 6 


-1.33526 


2.32719 


V S 


0.235331 


-4.21557 


Vio 


-0.025807 


-0.50636 



The corresponding zeros xo of this function, corresponding to the singularity of V(x), are 



L 


QM 


memb. 


1 


0.450158 


0.450158 


2 


0.506893 


0.506893 


3 


0.499717 


0.523646 


4 


0.500008 


0.514714 


5 


0.500000 


0.514469 



(40) 
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The value of a in each case is given by (2x ) 2 times the exact value (33) of a for QM, 



L 


QM 


memb. 


1 


0.0625 


0.0625 


2 


0.0792468 


0.0792468 


3 


0.0770189 


0.0845718 


4 


0.0771087 


0.0817113 


5 


0.0771062 


0.0816335 



(41) 



While the correct QM value (33) is approached very quickly, the convergence in the membrane case is slower. The 
fact that the last two values have such a small difference appears to be accidental. However, the results point towards 
a value above 0.080. 



VII. SUMMARY AND DISCUSSION 

We have used three methods to extract the pressure exerted by a tensionless membrane on two infinitely extended 
parallel walls from a five-loop calculation for smooth potentials. While variational perturbation theory with self- 
consistently determined q is converging too slowly for quantitative statements at five loops, variational perturbation 
theory with the assumption q = 1 gives a result a w 0.0813. The successive a values at the various loop orders in 
(24) suggest an extrapolated value of a between 0.0820 and 0.0825. Fixing a(g) to resemble the ^-structure of the 
ground state energy of a solvable quantum mechanics problem and analyzing the location of the singularities next to 
the origin of the resulting potential leads to a w 0.0816. As opposed to variational perturbation theory with q = 1, 
the sequence given by the considered loop orders gives only a modest indication of where a might settle. The results 
of both analyses point towards a value above the Monte Carlo result «mc = 0.0798 ± 0.0003 [5]. 

We have also studied the quantum mechanics problem with a potential for which we do not know the exact solution 
but which is very close to the potential of the solvable problem in the region of interest. We have investigated 
both variational perturbation theory with self-consistcntly determined q and with q = 1. The result is that the 
exponentially fast convergence of the solvable model towards the exact result for a cannot be expected for the general 
case, although good estimates of the exact result arc obtained. Since numerically, this case is close to what happens 
in the membrane case, this gives us an indication on how trustworthy our results are. As already noted at the end 
of sec. V, we can roughly state that the results for a in the quantum mechanics problem do not improve after the 
cll start to significantly increase in magnitude. Along this reasoning, one may still expect improving the results for 
the membrane problem by proceeding to higher loop orders. In particular, going to six loops appears feasible with 
reasonable effort, since only 8 of the 83 diagrams to be evaluated have no cutvertex and have therefore a true six- loop 
topology, as noted in Table II. 
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APPENDIX A: RECURSION RELATION FOR THE LOOP EXPANSION 

Here we define a loop expansion of the free energy and derive recursion relations for obtaining the required diagrams 
in a systematic way along the lines of [15]. We continue to work with k = ksT = 1. 
We write the energy functional (7) as 

E[p, G, {L^}} = -Gr 2 Vi¥>2 + L<°> + J2 L i 2fe) ,2^i ■ • • <P*k (Al) 

k=2 

with totally symmetric tensors G _1 and L^ 2k \ Their indices 1,2, ...,2k are shorthands for space arguments 
x\, . . . , X2k and a generalized Einstein convention implies integration over space arguments that appear twice in 
a term. Comparison with (7) shows that 

G^ 1 = 6 12 [dfdi + m 4 } (A2) 
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and 

4!. fe U = m 4 e2fcd 2(1 - fc) «i,...,2fc (A3) 

with 



d 2 xS(x — Xi) ■ ■ ■ S(x — x 2 k)- (A4) 



Note that the index of e 2 fc does not indicate a space argument and is exempted from the summation convention. 
The free energy Af m = — W is given by 

exp(^[G, {L^}}) = J DtpeM-Eif, G, {L^}}). (A5) 

W obeys the functional differential equation 

= J Dv^ivoexpi-EfaG^LW}])} 



v dG 02 <5G 02 dG 34 fe=3 &Ga2 S H,...,2kJ 



(A6) 



Splitting = T4 / | i( 2 fc)=0 + W/ = W + W7, so that VK obeys 

^01 + 2G ^g^=l) cx P(^o[G]) = 0, (A7) 

we get from (A6) 

where we additionally have identified indices and 1 and integrated over the respective variable. With 



and 



Eq. (A8) may be transformed into 



8G^ 2 



5 2 Wn 1 

(G 13 G 2i + G u G 23 ), (A10) 



<J12 J = -6LJ234G12G34 - 241,^234^12^35^46 777-^- ~ 8L1234G15G26G37G48 (77; — 77^ 1" 

OU12 OG56 \OG56OG78 OG56 OG78 



fc=3 

Define a loop expansion 



oc 

(2k) 

2k 

fe=3 



^1277751-^7 + ^11^22 



r r (2k-2) -11^22 I ( 2 fc-2) „ (2fe-2) 

,W -3 2,v \ ,)(j 12<'^3,...,2fe 12 dL 3,...,2k . 



(All) 



147 = 147(10 ( A12 ) 



L=0 



and set 



^ (0) - (A13) 
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by an appropriate normalization of the path integral measure Dp. Then 

W = W (1) = ~ HG- 1 )!! (A14) 

and 

oo 

W/ = E w(L) - (A15) 



L=2 



Eq. (All) separates into the two-loop equation 



SG 

and the recursion relation 



Gi2^7^ + 6L[% 4 G 12 G 34 = (A16) 



= 2 



78 



1.1- (2k) 
2 , KL, l,...,2k 

k=3 



G 12 (2fc-2) +2^11^22 (2fc _ 2) + ^ SG , (2k-2) 

dL 3,...,2k \ dLr 12 dL 3,...,2k 1=2 12 dL 3,...,2k , 



(A17) 



which holds for L > 2 and where we have taken into account that a diagram containing the tensor £( 2fe ) has at least 
k loops. 

Eq. (A16) is solved by 

W<® = -3L[% 4 G 12 G 34 . (A18) 

Before carrying out the recursion relation (A17), let us introduce a graphical representation of the resulting terms. 
Represent each free propagator G\ 2 by a line with two ends corresponding to the two space arguments x\ and x 2 and 
each tensor — n by a dot. Each line end is connected to the dot with an identical space argument. Then a dot 

corresponding to a tensor n has n line ends connected to it. In this way all terms appearing in the with 
L^l can be graphically represented. The zero-loop order is represented by a dot without lines, 

W {a) = • . (A19) 

Only does not fit into the graphical scheme above and as usual we use the graphical representation 

= lQ (A20) 

for it. Now we may write (A18) as 

^ (2)=3 QO< (A21) 

which is the starting point for the recursive determination of the other W^ L \ For instance, the three-loop contribution 
to W is 



w^ = 12 ^E) +36 GQO +15 rYy (A22) 



In Table II we list the numbers of different diagrams through seven loops and the numbers of diagrams at each loop 
order which have no cutvertex (by definition, upon cutting through such a vertex appropriately, a diagram decomposes 
into two diagrams, which consequently have independent momentum integrations) and have therefore their full loop 
topology. The contributions through five loops for both the QM and the membrane problem are collected in Table 
III. 
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number of loops L 


1 


2 


3 


4 


5 


6 


7 


diagrams 


1 


1 


3 


7 


24 


83 


376 


diagrams with L-loop topology 


1 





1 


1 


5 


8 


37 



TABLE II: Numbers of vacuum diagrams for some low loop orders and numbers of those with full loop topology. 



APPENDIX B: EVALUATION OF QUANTUM MECHANIC INTEGRALS 



All integrals in the QM case can be evaluated analytically. The propagator reads 

1 



A(k 2 ,m 2 ) 



k 2 + to 2 ' 



its Fourier transform is 



A(t,m) = 



dk e 



iki 



2tt k 2 + m 2 2m 



(Bl) 



(B2) 



For most integrals it is convenient to work in i-space, omitting the last ^-integration (which, due to time translation 
invariance, gives a factor T, the total interaction time). An exception are the one-loop integrals, which are computed 
easiest in momentum space. Using dimensional regularization, we get 



d D p d , , o 2\ 
D=1 DJ (2^V n(P +m) 



D=l 



d D p 



DJ (2n) D p 2 +m 



= to (B3) 



D = l 



for the only diverging integral, while the other one-loop integrals are given by 

f +oc dp 1 m^rtn-i) 

'J 71 



J — I 



2tt (p 2 + to 2 )" 2 ^T{n) ' 

Note the similarity to the membrane one-loop integrals (CI) and (C2) below, so that 

4 qm = Am- 1 J™", JT = 4m 2 "- 1 Jr mb - 



(B4) 



(B5) 



APPENDIX C: EVALUATION OF MEMBRANE INTEGRALS 

In the following, we will always assume that to 2 > 0. We give results in a form suited for numerical integration of 
the remaining loop momenta. 



1. One-Loop Vacuum Integrals 



The ubiquitous one-loop integrals without external momenta can be computed analytically as 



4 f^Wlx 41 1 [ d D P d 4 4 4 f d D p 



and 



P 

" 4 + TO 4 



D=2 



TO. 



2-4n 



r(n 



l p (p 4 + to 4 )" " 8 v^r(n) ' 
where dimensional regularization has been employed for J' . 



(C2) 
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2. One-Loop Bubble 

Several diagrams contain the one-loop bubble 

yQC = f A((k+p) 2 ,m 2 )A(p 2 ,m 2 ) 



4to 4 J p \ (fc + p) 2 + im 2 (fc + p) 2 — im 2 J \p 2 + im 2 p 2 — im 2 



4to 4 



with V'oi ' 1 ^ defined by 



Elementary integration gives 



(fc 2 , m?, m 2 .) _ ^ [(p + fc)2 + .^ ?](p2 + . mi) • 



Vfc 2 ± Aim 2 + Vk 2 
27r%/P'Vfc 2 ± 4im 2 Vfc 2 ± 4«m 2 - \/fc 2 



and 



so that 



^ 1) (fc 2 ,±m 2 , T m 2 ) = 



fc 2 + Vfc 4 - 4to 4 
47iVfc 4 - 4m 4 k 2 - Vk 4 - 4m 4 : 



dn 



1 



1 , k 2 + Vk 4 - Am 4 

In ■ — 2Re 



Vk 4 - Am 4 k 2 - Vk 4 - 4to 4 
Now we can also easily compute 



1 VA: 2 - 4im 2 + Vfc 2 

k 2 \/k 2 — Aim 2 \/k 2 — Aim 2 — Vk 2 



^ 2 \k 2 ,m 2 )= = f p A((k +P ) 2 ,m 2 )A(p 2 ,m 2 ) 2 = -^^A^\k 2 ,m 2 ) 



5fc 2 



47rm 4 (fc 4 - 4m 4 )(/c 4 + 16m 4 ) 
1 



(C3) 



(C4) 



(C5) 



(C6) 



(C7) 



167rm 8 



k 4 - 6m 4 , fc 2 + Vfc 4 -4m 4 

In , — 2Rc 



(fc 4 -4m 4 ) 3 / 2 fc 2 - Vfc 4 -4m 4 



fc 2 + him 2 , \/fc 2 +4jm 2 + 
In 



fc 2 (fc 2 +4«TO 2 ) 3 / 2 V^ 2 +4im 2 - 



fc 2 
P 



(C8) 



3. Sunset Self-Energy 

For several diagrams we need to compute 
A ss (fc 2 ,TO 2 ) = = f A((k+p) 2 ,m 2 )A((p + q) 2 ,m 2 )A(q 2 ,m 2 ) 

Jpq 

= ~^—e [VWfc 2 ,m 2 ,m 2 ,to 2 ) - ip ss {k 2 , m 2 , to 2 , -to 2 ) - ip ss (k 2 , to 2 , -m 2 , to 2 ) + ip ss (k 2 ,m 2 , -m 2 , -to 2 ) 

— lp ss (k 2 , —TO 2 , TO 2 , TO 2 ) + 1p ss {k 2 , —TO 2 , TO 2 , —TO 2 ) + 1p ss {k 2 , —TO 2 , — TO 2 , TO 2 ) — l/> ss (fc 2 , —TO 2 , —TO 2 , —TO 2 )] 

(C9) 



with V'ss defined by 



V> ss (A: 2 ,to 2 ,TO2,TO3) = / 

Jpq 



[(fc + p) 2 + zto 2 ] [(p + g) 2 + im 2 ] (q 2 + ito§) ' 



(CIO) 
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Note that 



ip ss (k 2 , m 2 , to 2 , m 2 ) = ip ss (k z , —mi 2 , —to 2 , —to 



2\* 



and 



^(fc 2 , -to 2 , to 2 , to 2 ) = t[; ss (k 2 ,m 2 ,~m 2 ,m 2 ) = ip ss (k 2 , to 2 , to 2 , -to 2 ) 
= tp ss (k 2 , to 2 , —to 2 , —to 2 )* = ip ss (k 2 , —m 2 ,m 2 , —to 2 )* = tp ss (k 2 , —to 2 , —to 2 , to 2 )*. 



(Cll) 



(C12) 



ip ss may be evaluated as 

i/> ss (fc 2 ,TO 2 ,TO2,TO 2 ) = da 

Jo Jp 



(k + p) 2 + im\ J„ [q 2 + 2ap ■ q + ap 2 + aim 2 , + (1 — a)im 



212 



^ la 



[(k +p) 2 + zto 2 ][q:(1 - a)p 2 + aim 2 , + (1 — a)im\] 



da 



4tt J a a(l - a) J 



dp 



P 2 + 2(1 - f3)p ■ k(l - p)(k 2 + im 2 ) + pi(^ + 



= {**?L a(l-a)J ^ m _ 



13(1 - P)k 2 + (1 - /3)zto 2 + (3i + ^) 



da / d/3 



(4tt) 2 7 7o a(l-a)(l-/3)(/?A: 2 + iTO 2 )+/3i[aTO| + (l-a)TO 2 ] 



(4tt) 2 7 

If to 2 . = to|, this becomes 

tp ss (k 2 ,m\,ml,ml) -- 



+i 



d(3 j dx 



(1 - x 2 ){\ - /3)(/3fc 2 + zto 2 ) + 2i/3[(l + x)to 2 + (1 - x)to 2 



(C13) 



+i 



d/3 / c?x 



(4tt) 2 7 (1 -x 2 ){\ -(3)(l3k 2 +im 2 ) + Aipm 



2 " 



dp 



(4tt) 2 J (l-P)(pk 2 +im 2 ) J_, 



+i 



dx- 



1 + 



4i.dj, 



(C14) 



(l-/3)(/3fc 2 +imJ) 



while for to 2 = — to?, we have 



V> ss (fc ,TO 1 ,TO 2 ,-TO 2 ) 



+1 



d/3 / dx 



(4tt) 2 7 (1 -a; 2 )(l -/3)(/3fc 2 +im 2 ) +4ia;/3TOi 



2 ' 



(4tt) 2 7 (l-/3)(/3fc 2 +im?) 7-i i _ x 2 + 



+i 



4if3m* 



(C15) 



(l-/3)(/3fc 2 +m^) 



Noting that the formulas 



h{z)= L v^ = ^L dx hh-x + 7^ x ) 



ln(yi+l)-ln(VI-l) 



(C16) 



and 

,+i 

AW -L r 

1 



dx 



x 2 + 2zx 2\/l 



1 f +1 , ( 1 
/ dx ■ h 

+ z 2 J-i + z 2 - z + x 



Vl + z 2 + z - . 



ln( VTTI 2 - z + 1) - ln( VTTI 2 - z - 1) - ln(\/i + z 2 + z - 1) + ln( VT+1 2 + z + 1)1 (C17) 



are numerically safe to use if the branch cut of the logarithm is taken from to —00, we may now express the V'ss in 
(C9) with the help of /3-integrals involving /1 and fi- However, these integrals are difficult to evaluate numerically if 
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m l = m i = — m i- These cases may be avoided by making use of (C12). Together with (Cll) we get 



A ss (fc 2 , to 2 ) = ^-glm[2^ ss (fc 2 , to 2 , m 2 , to 2 ) - 6ip ss (k 2 , to 2 , to 2 , -to 2 )] 



dp 



32tt 2 to 6 J (l-f3)([3k 2 +im 2 ) 
This form is easy to implement numerically. 



h i + 



4?/3to 2 



(l-/3)(/3A: 2 + iTTi 2 



3/ 2 



2i(5m 2 



(l-f3)(f3k 2 +im 2 ) 



(C18) 



4. Eye Bubble and Triangle Coupling 

For the integral 15-3, we need the eye bubble subdiagram 

A cyc (fc 2 ,TO 2 )= 3(0); =^A tr (fc,r,m 2 ) 2 , 
where the triangle subdiagram A tr is given by 



A tr (fc,r, to 2 ) 



-p+f*^ ~ J p 



"8(27t) 2 to6 



1 



(4tt) 2 to 6 



with xpt r defined by 



A(p 2 , m 2 )A((p + k) 2 , m 2 )A((p + r) 2 , m 2 ) 
\iptr{k, r, to 2 , to 2 , to 2 ) — tptr(k, r, to 2 , to 2 , —to 2 ) 

— 1ptr(k, T, TO 2 , —TO 2 , TO 2 ) + tptr(k, T, TO 2 , —TO 2 , —TO 2 ) 

— tptr(k, r, —m 2 ,m 2 , to 2 ) + iptr{k, r, —to 2 , to 2 , —to 2 ) 

+ 1ptr(k, T, — TO 2 , —TO 2 , TO 2 ) — 1p\, v {k, T, — m 2 , —TO 2 , —TO 2 )] 

Im [tpt T (k, r, to 2 , to 2 , to 2 ) — tptr(k, r, to 2 , to 2 , — m 2 ) 

— ipti{k, r, to 2 , —to 2 , to 2 ) + iptr(k, r, to 2 , —to 2 , —to 2 )] 



(i 2 p 



(p 2 + im 2 ) [(p + fc) 2 + im 2 ] [(p + r) 2 + im 2 ] 

1 



/ dpp I in /- — -j - — - - - — — — - --— - - - — - - — 

J J q (p 2 + im\)[p 2 + k 2 + 2pkcos((f> - fa) + im\]\p 2 + r 2 + 2pr cos(0 - cf> r ) + im 2 ] 



/•2-7T pC 

/ # / 

Jo Jo 



dpp 



(P - a+)(p - a-)(p - b+)(p - 6_)(p - c+)(p - c_) 



1 



/o a + -a_ 



a_L In Oj 



a_ In a_ 



+ 



1 



+ ■ 



1 

c + — C- 



(a + — 6 + )(a + — 6_)(a + — c + )(a + — c_) (a_— 6 + )(a_— 6_)(a_— c + )(a_— c_) 

fe + ln& + fo_ln6_ 

(& + -a + )(&+-a_)(& + -c + )(& + -c_) (6_-a+)(&_-a_)(6_-c + )(6_-c_) 

c + In c + c_ In c_ 



(c+-a+)(c+-a_)(c + -6 + )(c + -6_) (c_-a + )(c_-a_)(c_-& + )(c_-&_) 



with 



a± = ±iy imf 



(C19) 



(C20) 



(C21) 



(C22) 



b± = — k cos((/) — <f>k) ± iy k 2 sin {<p — 4>k) + irn^, 



(C23) 
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c± = —rcos(4> — <f> r ) ± i\jr 2 sin 2 (0 — <j> r ) + im 2 . (C24) 
For the numerical evaluation of ^5-3, it is also useful to know the large-fc behavior 

A eyc (fc 2 ,m 2 ) - A /3 _ l7 fc 2 »m 2 (C25) 

with from Table III. 

Let us remark here that V'tr and therefore A tr can also be computed analytically (e.g. by performing the integrals 
over the cartesian components p x and p y of p), but the resulting expressions are rather lengthy and it may be a 
delicate issue to remain on the same Riemann sheet while evaluating the logarithms and square roots involved. 

5. Numerical Considerations 

Most integrals are evaluated rather easily. They are either known analytically involve only one numeric integration 
or involve one such integration involving the function A ss (fc 2 ,m 2 ), whose evaluation implies one numeric integration, 
see (C18). These cases are easily dealt with by using any standard software integration package, e.g. mathematica, 
which was used here. A lot of integrals can be evaluated in different ways, so the safety of using A ss (fc 2 ,m 2 ) 
inside an integral can be and has been checked. For instance, the integral I3-1 can be performed by integrating 
either A^}' 1 ^ (k 2 , m 2 ) 2 or A ss (fc 2 , m 2 )A(k 2 , m 2 ) over k 2 . By such cross checks and by varying the settings within 
mathematica for the numerical integrations, one may easily gain confidence that the integrations performed are 
accurate through the number of digits given in Table III. 

The only exception to the above considerations through five loops is integral 15-3. As indicated in Table III, the 
product A^j' 1 ^ (fc 2 , m 2 )A oyc (fc 2 , m 2 ) has to be integrated over k 2 . While A^}' 1 ^ {k 2 , to 2 ) is known analytically from (C7), 
A e y e (fc 2 ,TO 2 ) involves a further two-dimensional numeric integration of A tr (fc,r, to 2 ) 2 over r, see (C19). A tr (fc,r, to 2 ) 
itself implies a one-dimensional numeric integration, see (C20) and (C21). Using this method, an uncertainty in the 
last digit given in Table III remained in our computations. We achieved the precision through the last digit given in 
Table III by rewriting / 5 _ 3 as the five-dimensional integral 

1 n Jo Jo Jo Jo Jo 

x A(p 2 , to 2 )A(<z 2 , m 2 )A((p + r) 2 , m 2 )A((q + r) 2 , m 2 ) A^ (r 2 , to 2 )A^ 1} ((p - q) 2 , to 2 ) (C26) 

in an obvious notation and asking mathematica for successively increasing precisions of the final result. 

The precision achieved for all integrals is several orders of magnitude better than needed to determine, say, the 
first three non-zero digits of our estimates for a. Higher precision is easily attained for all integrals but I5-3, but is 
unnecessary for the purpose of this work. 



17 



TABLE III: Diagrams L—n (nth L-loop diagram) through five loops 
and their combinatorial factors cl-u, coupling constant factors gh-n 
and values Il-u of the corresponding integrals for m = 1. D — 1 
corresponds to the QM problem and D — 2 to the membrane problem. 



0-1 



L—n 


diagram 


CL-n 


Qh-n 


I L - n for D = 1 


7 L _ n for D = 2 



-£o 



1-1 



o 



1/2 



-Jo = -1 



-Jo = -1/4 



2-1 



oo 



-£4 



1/4 



J 2 = 1/64 



3-1 



3-2 




QOO 



12 



36 



4 



e\ 



1/32 
1/16 



/ t Ai; 4) (t> 2 f = 4.04576 x 10" 
Ji 2 J 2 = 1/1024 



3-3 




15 



1/8 



J't = 1/512 



4-1 

4-2 
4-3 

4-4 





OQQO 




288 

576 
432 

288 



-el 



-el 



-el 



-el 



3/512 

5/512 
1/64 

3/128 



J k A^' 1 ^ (fe 2 ,m 2 ) 3 = 1.63237 x 10" 

fJi/s-i = 3.16075 x 10" 
jfjf = 1/16384 

Jfj 3 = 3/32768 



4-5 



4-6 





360 



540 



£4 £6 



£4 £6 



1/64 
1/32 



Jih-i = 5.05719 x 10" 



J1J2 = 1/8192 



4-7 




105 



-£8 



1/16 



Jf = 1/4096 



5-1 



5-2 



5-3 



o 




2592 



2304 



10368 



£4 



4 

£4 



e\ 



5/4096 



19/12288 



7/6144 



/ A^ik 2 ,™ 2 ) 4 = 7.55133 x 10- 7 



/ A(fc 2 ,m 2 ) 2 A as (fc 2 ,m 2 ) 2 = 1.04187 x 10" 



J k A^' 1) (fc 2 ,m 2 )A cye (fc 2 ,m 2 ) = 6.71540 x 10" 7 
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TABLE III: (continued) 



L—n 


diagram 


CL-n 


9L-n 


I L - n for D = 1 


I L _ n for D = 2 



5-4 



5-5 





5-6 




5-7 



5-8 



5-9 



5-10 




OQ3GO 





20736 



10368 



6912 



6912 



5184 



10368 



2592 



§Ji/ 4 -i = 1.36031 x 10~ 6 



J 2 f fc A^ 2) (fc 2 ,m 2 ) 2 = 2.54723 x 10" 



J\ J k A(fc 2 ,m 2 ) 3 A ss (fc 2 ,m 2 ) = 3.09329 x 10" 



f J1J2/3-1 = 1.97547 x 10~ 6 
J1J2 = 1/262144 

Ji J2J3 = 3/524288 
jfJi = 5/524288 



5-11 



5-12 



5-13 



5-14 



5-15 



5-16 









5760 



12960 



4320 



17280 



4320 



6480 



7/3072 



3/1024 



5/1024 



J fc A(fc 2 ,m 2 )A ss (fc 2 



3/256 



J1J4-1 = 2.04047 x 10~ 6 
|Ji 2 /3-i = 3.95093 x 10~ 6 

|Ji 2 /3-i = 3.95093 x 10~ 6 
J1J2I3-1 = 3.16075 x 10~ 6 

jfJz = 3/262144 
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TABLE III: (continued) 



L—n 


diagram 


CL-n 


9L-n 


I L - n for D = 1 


I L _ n for D = 2 



5-17 



5-18 





6480 



6480 



1/128 



1/128 



J3j| = 1/131072 



J3j| = 1/131072 



5-19 



5-20 



5-21 






360 



2700 



2025 



4 



4 



1/192 

1/128 
1/64 



/ A ss (fc 2 ,m 2 ) 2 = 3.76084 x 10~ 



Jth-i = 6.32149 x 10~ 



Jf J 2 = 1/65536 



5-22 



5-23 





5040 



5040 



£4^8 



£4^8 



1/128 



1/64 



Jlh-i = 6.32149 x 10~ 6 



Ji 4 J 2 = 1/65536 



5-24 




945 



-eio 



1/32 



j\ = 1/32768 
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